Single-channel transmission in gold one-atom contacts and chains 
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We induce superconductivity by proximity effect in thin layers of gold and study the number 
of conduction channels which contribute to the current in one-atom contacts and atomic wires. 
The atomic contacts and wires are fabricated with a Scanning Tunneling Microscope. The set of 
transmission probabilities of the conduction channels is obtained from the analysis of the /(V) 
characteristic curve which is highly non-linear due to multiple Andreev reflections. In agreement 
with theoretical calculations we find that there is only one channel which is almost completely open. 
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Much effort has been devoted in the last decade to the 
understanding of electron transport processes and me- 
chanical properties of atomic-sized point contacts and 
chains between metallic electrodesi. Coherent electron 
transport in these nanostructures can be understood in 
the frame of the scattering formalism. The conduc- 
tance G of these nanocontacts is given by the Landauer 
formula G = Go X^ti T ii where Tj are the transmission 
probabilities for each of N conductance channels and 
Go = 2e 2 /h is the conductance quantum. For a given 
contact realization the conductance channels are in gen- 
eral neither completely open nor completely closed and 
the transparency Tj of each channel depends on the mate- 
rial forming the contact 2 ^, detailed atomic arrangement 
and applied stress^. In this report we study the chan- 
nel transparency set {t^} both in gold one-atom contacts 
and chains of single gold atoms 5 , for which theoretical 
models^iLSiS p re dict an almost completely open single 
channel, and therefore a conductance close to Go- 

It is well established that conductance histograms 
(CH) show that the conductance of one-atom contacts 
of gold is close to Gcr 1 but the values of Tj cannot be ob- 
tained from the only measurement of the total conduc- 
tance. However, the marked non-linearity of the current- 
voltage characteristic (IV) of superconducting contacts 
has been exploited to obtain the set of transmission co- 
efficients {ri} of atomic-sized aluminum contacts^. The 
channel transmission probabilities are extracted by fit- 
ting the measured IV curve to a sum of N independent 
IV s calculated for individual channels with a given trans- 
mission probabilitjiiiiiS. This method was afterwards ex- 
tended to other superconducting materials, showing that 
the number of conducting channels contributing to the 
current in one-atom contacts is limited by the valence of 
the atom at the contact. 

Gold is not superconducting. However we take ad- 
vantage of the superconducting correlations induced in a 



thin layer of normal metal in contact with a supercon- 
ductor. The energy spectrum is modified and a gap is 
opened at the Fermi energy. At low voltages transport 
is dominated by Andreev reflection processes which re- 
sults in non-linear IV curves. Here we use the aforemen- 
tioned method 2 *^ to analyze the transmission coefficients 
of proximity induced superconducting one-atom contacts 
and atomic chains of gold. 

The proximity effect has been previously exploited 2 *^ 
to get information of the {Tj} in gold atomic contacts. 
In the first experiment the IV curves were fitted to 
a sum of theoretical IVs corresponding to BCS super- 
conductors and the contribution of a single channel in 
one-atom contacts is reported. However, the energy de- 
pendence of the density of states and the probability of 
Andreev reflection at a proximity induced superconduct- 
ing contact are smeared with respect to the BCS model 15 . 
The modifications due to the proximity character of the 
superconducting correlations were taken into account in 
the later experiment^ 3 .. Most of the experimental IVs 
recorded at the last conductance plateau before contact 
breaking could be fitted with a single channel. Several 
channels were necessary to fit the IV curves of some con- 
tacts with conductance smaller than Go- However, in 
these experiments the conductance of the smallest con- 
tacts is usually much smaller than Go, being the first 
peak of the histogram placed at 0.6 Go and plateaux not 
nearly flat. These results are not consistent with what 
has been usually reported in gold atomic contacts in the 
normal state*. Here we report experiments in which the 
conductance histograms for proximity superconducting 
gold and normal gold are in remarkable agreement (Fig 
1). We also study the channel content of chains of gold 
atoms. 

Nanocontacts are produced by pressing two wires 
crosswise against each other. The wires are used as elec- 
trodes in place of tip and sample of an STM. The separa- 
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FIG. 1: Differential conductance at a bias voltage of f0 mV 
of proximity induced superconducting gold contact during 
pull-off. One-atom contact formation (left) and atomic chain 
(right). The length of the conductance plateau indicates that 
a chain longer than five atoms is formed. Inset: conductance 
histograms of gold in the normal state (bottom) and proxim- 
ity superconducting (shifted upwards 10000 counts). 



tion and contact size between the wires can be controlled 
with the piezoelectric positioning system the microscope. 
The advantage of using two crossed wires is the possibil- 
ity of selecting the position of the point contact along 
both wires using the coarse lateral displacement capa- 
bility of the microscope, allowing for the exploration of 
point contacts at different spots along the wires. The 
wires (0.25 mm diameter) are made of bulk lead and are 
in a first preparation step covered by thermal evaporation 
with a thick layer of lead (900 nm at a rate of 0.8 nm/s). 
This thick layer of lead provides a clean surface. A thin 
layer of gold (28 nm, rate of 0.1 nm/s) is then evaporated 
on top of the Pb layer. The lead and gold deposition se- 
quence is performed without breaking the vacuum in the 
chamber (< 10~ 6 mbar), thus preventing the presence of 
oxide at the Pb-Au interface and ensuring good electri- 
cal contact between both layers. The substrate is at a 
temperature of 80 K during film deposition. One-atom 
contacts are fabricated by slight indentation (< 3 nm) 
and subsequent retraction of the electrodes. Further re- 
traction results in contact breaking and a jump to the 
tunnelling regime. Experiments are done at 1.8 K. Dur- 
ing nanocontact pull-off the conductance trace is step- 
like (see Fig. 1). This characteristic behavior has been 
shown to be due to the mechanical processes that take 
place during contact breaking 16 ; conductance plateaus 
corresponding to elastic deformation stages and sharp 
conductance changes related to sudden rearrangements 
of the atoms in the narrowest part of the nanocontact. 

For gold nanocontacts in the normal state the last con- 
ductance plateau before contact breaking is close to Go 
and corresponds to the smallest possible contact: a one- 
atom contact. Despite the inherent variability of the ex- 



act conductance trace during contact pull-off, it has been 
shown in several experiments that the conductance his- 
togram (CH) over a large number of contact breaking re- 
alizations displays peaks at conductance values that are 
characteristic of the chemical nature of electrodes^. CH 
in gold nanocontacts have a characteristic first peak at 
a conductance close to Go- We show in the inset of Fig. 
1 a comparison between CH obtained for bulk gold tips 
in the normal state and in proximity induced supercon- 
ducting gold nanocontacts. Due to the presence of excess 
current in the IVs of superconducting nanocontacts it is 
necessary to measure the differential conductance at a 
fixed bias voltage well above A/e, where A is the energy 
of the superconducting gap. The remarkable agreement 
between both CH supports the validity of our analysis 
also for atomic-sized contacts of gold in the normal state. 
The low value of the conductance of the first peak in the 
conductance histogram in Ref. was probably due to 
enhanced elastic scattering related to sample preparation 
method. 

We show in Fig. 2 representative IV curves recorded 
at a one-atom contact, in an atomic chain and in the tun- 
neling regime, up to a voltage 2Ao/e, where Ao = 1.40 
meV is the bulk superconducting gap for lead. We show 
both the measured TV curves (symbols) and theoretical 
fitting (lines). The theoretical IVs in the contact regime 
are calculated by solving the time-dependent Bogoliubov 
de Gennes equations within the scattering formalismii*i&. 
This method requires the knowledge of the Andreev re- 
flection amplitude of probability (a(E)) at the contact, 
where the voltage drops. a(E) can be computed if the 
normal and anomalous Green functions are known— In 
the tunneling regime the IVs are calculated by the usual 
convolution of the densities of states at both sides of the 
barrier—. 

As both Au-Pb electrodes are only very weakly coupled 
at the weak link, to calculate the density of states and 
the Andreev reflection amplitude of probability, we model 
our system as two independent normal-superconducting 
(NS) structures and solve self-consistently the Usadel 
equation a 14 ! 15 . Usadel equations provide a quasiclassical 
description of the Green functions of a superconductor in 
the dirty limit, in which electronic transport is diffussive. 
Elastic impurity scattering is included in the Born ap- 
proximation and is characterized by a mean free path 
I (or a diffusion coefficient D — vpl/3). A description 
based on Usadel equations was recently used to explain 
the IV curves of lead nanostructures under the influence 
of a magnetic field and proximity effect, providing excel- 
lent quantitative agreement with experiment both in the 
contact and tunneling regim o 19 ! 20 . 

The two NS are assumed to be equal and consist of a 
dirty normal layer with thickness d,N which is bounded 
at one end by vacuum and joined to a dirty semi-infinite 
superconductor at the other. Note that the diffussive 
description is expected to be valid in the case "C I- 
The Green functions which are relevant to calculate the 
IV curves are the ones at the vacuum-bounded edge of 
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FIG. 2: Main figures show the IV curves corresponding to the 
tunneling regime (b), a single-atom contact (bottom curve in 
(a)) and an atomic chain (top curve in (b)). The curve corre- 
sponding to the atomic chain has been displaced for clarity. 
Experimental results are shown by symbols and theoretical fit- 
tings by solid lines. The transmissions obtained in the contact 
regime are T\ = 0.995 (single atom), and T = 0.96 (chain). 
Inset in (a) shows the density of states used to calculate the 
theoretical curves -see text. Inset in (b) show the derivative 
of the experimental (the one with noise) and calculated tun- 
neling curves. 



the normal layer. The proximity effect is also, affected 
by the existence of a barrier at the interface, as quasi- 
particles normally reflected do not contribute to it. In a 
diffusive system, the mismatch between the characteris- 
tic parameters (conductivity and diffusion coefficients) of 
the normal and superconducting metals, leads to an ef- 
fective barrier for the quasiparticles. It can be described 
through the parameter V = (Dg/ Djv) 1/,2 <7jv /&s- In our 
model we assume T = 1 and a vanishing resistance of 
the interfaceSi. With these assumptions the supercon- 
ducting correlations in the normal metal are described 
by A and the value of djv/£s, where £s is the coher- 
ence length of the superconductor. In the atomic chain 
and one-atom contact the transmission channels set en- 
ters also as fitting parameters. A similar model was used 
by Scheer et in the description of Au contacts with 
superconductiviy induced by proximity effect. In their 
case, however, the mismatch parameter T was also used 
as a fitting parameter. 

Solid lines in Fig. 2 show the calculated IV which 
for the same sample fitting parameter, (In/£,s = 0.81, 
best fit both the curves in the tunneling and the contact 
regime. Assuming the nominal thickness = 28 nm it 
corresponds to a superconducting coherence length £s — 
34 n m, in good agreement with the values obtained in 
Ref. GJ and |M 

The tunneling regime IV curve, plotted in Fig. 2b, 
shows a gap smaller than Ao and a bump characteristic 
of NS structures 2 ^. It is however poorly fitted^, which 



means that the density of states used, shown in the inset 
of Fig. 2(a), does not agree completely with the experi- 
mental one. The source of disagreement in the fitting can 
be better understood by looking at the derivative of the 
experimental and theoretical tunnel IV curves, shown 
in inset in Fig. 2 (b). The theoretical curve shows an 
asymmetric peak at the gap which decays slowly as the 
energy increases and results from the one in the density 
of states. This asymmetric peak structure is character- 
istic of the diffusive regime for parameters which give 
gaps in the normal metal much smaller than Ao- The 
peak in the experimental curve is more symmetric. More 
rounded peaks (and tunneling curves with a bump more 
similar to the one obtained experimentally) can be ob- 
tained for much smaller values of d^/^s- I n the diffusive 
regime the gap induced in the spectrum of the normal 
metal decreases with djv/Cs- Lower values of the fitting 
parameter induce gaps with values much closer to Ao and 
could not explain the strong reduction of the gap found 
experimentally and would give worse fits. 

Scheer et alm& also found disagreement in the fit of the 
tunneling conductance and related it to the non-diffusive 
character of the samples, being the condition djy <SC I not 
fulfilled. The proximity effect is a consequence of the co- 
herent superposition of Andreev reflection and is strongly 
influenced by the length of the path which the electron 
travels between Andreev reflection processes. Thus, the 
energy dependence of the induced pair correlations is very 
sensitive to the degree of disorder in the normal metal 
and the shape of the spectrum differs considerably in the 
clean and dirty limitsi^SSiS. We have also tried, with- 
out success, to fit the experimental curves assuming that 
there is no disorder being the transport in the electrodes 
ballistic, instead of diffusivo 22 i 23 i 24 ' 25 . As concluded by 
Scheer et alm&, we think that the lack of good fittings in 
the tunneling regime is due to transport in the gold layer 
being in the weakly disordered regime, instead of in the 
diffusive or clean limits. 

Given the uncertainty in the density of states, the ac- 
curacy with which the channel transmission set is ob- 
tained in the contact regime is slightly reduced, com- 
pared to the ones done in BCS superconductors. The 
single-atom curve can be reasonably well fitted with a 
very open channel with transmission T± = 0.995. The 
chain is reasonably well fitted by a single channel with 
transmission T — 0.96. This result is representative of 
the general behavior that we observe thus confirming that 
one widely open channel is responsible for the conduction 
in single atom gold contacts, in agreement with theoret- 
ical predictionsS^SiS. 
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